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ABSTRACT 

The Tycho supernova remnant (SNR), as one of the few historical SNRs, 
has been widely studied in various wavebands. Observations show evidence that 
Tycho is expanding in a medium with density gradient and possibly interacting 
with dense ambient medium toward the northeast direction. From the FCRAO 
CO survey of the outer Galaxy, we have identified a patch of molecular clouds in 
this area and have conducted a follow-up observation with Nobeyama 45m radio 
telescope. The high-resolution (16") Nobeyama data shows that a large molecular 
cloud surrounds the SNR along the northeastern boundary. We suggest that the 
Tycho SNR and the molecular cloud are located in the Perseus arm and that 
the dense medium interacting with the SNR is possibly the molecular cloud. 
We also discuss the possible connection between the molecular cloud and the 
Balmer-dominated optical filaments, and suggest that the preshock gas may be 
accelerated within the cosmic ray and/or fast neutral precursor. 

Subject headings: supernova remnants - ISM: individual (G120. 1+1.4) - ISM: 
molecules 



1. INTRODUCTION 



The Tycho supernova remnant (SNR) is known as a remnant of a Type la supernova 
that occurred in the year 1572 (Baade 1945). As one of the few historical SNRs, the rem- 
nant has been widely studied in various wavebands. In radio continuum, it shows a shell 
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morphology with clearly defined outer boundary and clumpy inner structure (e.g., Dickel 

et al. 1991). Recent high resolution X-ray images also exhibit similar features (Decourchelle 
et al. 2001; Hwang et al. 2002). Chevalier et al. (1980) derived a shock velocity by modelling 
the optical spectrum from a fast nonradiative shock and compared it with an optical proper 
motion (Kamper & van den Bergh 1978), giving an estimated distance of 2.3 ± 0.5 kpc. 
Other independent estimates give similar distances (Albinson et al. 1986; Strom 1988). It 
is, however, noted that there has been controversy in interpreting the H I absorption data, 
giving a different distance (e.g., Schwarz et al. 1995). This will be discussed in § 3.1. 

Observations show evidence that Tycho is interacting with ambient dense clouds toward 
its northeast (NE) direction. Reynoso et al. (1997) has studied its expansion using a 10- year- 
lap VLA data and found that the NE area, which is farthest from the geometrical center of 
the SNR, is currently moving slowest, suggesting that the remnant is expanding into a higher 
density medium. Reynoso et al. (1999) has conducted a H I absorption study towards the 
remnant and suggested that Tycho is currently interacting with a dense H I concentration 
in the NE (wlsr = —51.5 km s^^ ) which locally slows down the expansion of the shock 
front. There are other observations that also suggest its interaction with inhomogeneous 
environment (Velazquez et al. 1998; Katz-Stone et al. 2000; Douvion et al. 2001). 

From the FCRAO ^^CO survey of the outer Galaxy (Heyer et al. 1998), we have identified 
molecular clouds which could be associated with Tycho (Figure 1). They are arc-shaped and 
the inner boundary roughly traces the NE rim of Tycho. The morphology of the molecular 
clouds suggests physical association with Tycho. We have conducted high-resolution ^^CO 
J=l — line observations and present the results in this Letter. 



2. OBSERVATION AND RESULT 

The observation was carried out for a total of 20 hours during 2003 January 11-13 using 
the Nobeyama 45m radio telescope. The BEARS focal-plane receiver system was used to 
cover 12' x 12' area centered at Tycho whose diameter is ~ 8'. HPBW of the 45m telescope at 
^^CO (J=l-0) is about 16" and we have mapped the region with a 8" grid spacing. Pointing 
accuracy was quite good during the observing run and we estimate it to be better than 
5". Velocity resolution of 0.08 km was achieved using 1024 channel autocorrelator as 
a backend. The data were calibrated using the standard procedure with NewStar software 
developed by NAOJ. Exposure time is not uniform, leading to varying noises. The NE 
portion of the image has the most significant S/N (ATj-ms ~ 0.1 K) while the inner part of 
the remnant has a poorer S/N ratio (ATrms — 0.3 K). 
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In Figure 2, we present channel maps in gray scale together with the radio continuum 
showing the boundary of Tycho. There is virtually no ^^CO emission at wlsr > ~55 km s~^ 
(except from the gas in the local arm) and most of the emission is between —67 km s~^ < 
■^LSR < ~60 km s~^ with some faint emission localized in the northern part of the remnant at 
vlsr = —58 ~ —55 km s^^. The emission is generally from regions surrounding the remnant. 
In particular, at velocities between —63.5 to —61.5 km s~^, the emitting area appears to be 
in contact with the remnant along its NE boundary. This is more clearly seen in Figure 3, 
the integrated intensity map over i>lsr = —63 ~ —60 km s~^. In particular, at the location 
of optical knot g where the radio continuum boundary is deformed and the radio expansion 
is slowest, molecular clumps are touching the radio boundary. 

3. DISCUSSION 

The ambient environment of Tycho is not homogeneous as clearly revealed by the slower 
expansion of the NE rim. The morphology of the molecular cloud suggests their physical 
interaction. We discuss their physical association in the following sections. 

3.1. Distance to the Tycho SNR and the Molecular Cloud 

The distance to Tycho has been estimated by several authors using independent methods 
(Kamper & van den Bergh 1978; Albinson et al. 1986; Strom 1988). Most estimates agree 
on a distance of ~ 2.3 kpc except for the estimate based on H I absorption observations. 
Applying the H I absorption technique to Tycho is not straightforward due to the non- 
circular motion of the gas in the Perseus arm. The peculiar velocity field of the Perseus arm 
is well known, e.g., radial velocities of H II regions and OB stars in the Perseus arm are 
lower than those derived from photomctric/spcctroscopic distances (sec Schwarz et al. 1995, 
and reference therein). The widely accepted model is the "two armed spiral shock" (TASS) 
model proposed by Roberts (1972). In this model, the Galactic spiral pattern is visualized 
by a large scale shock residing in the Galactic density wave. Multi-value velocity-distance 
relation is naturally expected in this scenario, and, toward Tycho {i — 120°), the radial 
velocity is expected to drop abruptly from —20 km s"^ to —50 (±8) km s~^ at the shock 
front at 2.2 kpc and slowly increase to —40 km s~^ at 2.8 kpc, after which it slowly decreases 
(see Figure 14 of Schwarz et al. 1995). 

The velocity of the most negative absorption feature towards Tycho is around —50 km s~^ 
with weak absorption at —59 km s~^. Albinson et al. (1986) attributed the absorption to 
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the gas behind the spiral shock and placed Tycho just behind the large scale spiral shock, 
at a distance of 2.21q5 kpc. They noted that there is little absorption between —45 and 
—38 km s~^ which was considered as evidence for Tycho not being far beyond the spiral 
shock front. On the other hand, Schwarz et al. (1995) has argued that it should be placed 
at a distance of 4.6 ± 0.5 kpc where almost circular rotation is applicable. One of their 
arguments was that the velocity of the most negative absorption feature (—59 km s~^) is too 
low compared to the model-predicted shock velocity (—50 km s~^). The velocity structure 
in the Perseus arm, however, is rather highly dependent on the Galactic latitude. This trend 
can be clearly seen from the recent FCRAO CO survey of the outer Galaxy (Heyer et al. 
1998) and DRAO low resolution H 1 survey (Higgs Sz Tapping 2000). The peak velocity 
of ^^CO associated with Perseus arm in the region of / ~ 120° is around ~ —50 km at 
b ~ 0°, but systematically decreases to ~ —60 km s~^ at 6 ~ 1° and further decreases for 
higher latitudes. If these ^^CO molecules arc associated with the high density regime behind 
the spiral shock front, as is likely to, the shock velocity itself should have a similar trend, 
although it is not explicitly accounted for in the TASS model. H I data can be described in 
a similar way. Therefore, if we account the 6-dependent behavior of the velocity field, the 
absorption velocity could be as low as —60 km s~^ even if Tycho is just behind the shock 
front, and we consider that the H I absorption experiments are consistent with Tycho being 
at the near side of the Perseus arm. 

The velocity of the CO molecular cloud in the NE area of Tycho is ~ —62 km/s, which is 
slightly lower than the minimum velocity of the H I absorption feature. The cloud, and other 
clouds in Figure 2, however, seem to be part of a large-scale distribution of molecular gas 
presumably representing the high-density ridge behind the spiral shock. The CO emission 
between £ = 119° and 121° is generally confined within a velocity interval of —66 and 
—60 km s^^. This is much narrower than that of the H I gas (—70 ~ —45km s^^), whose 
peak velocity is also blue-shifted from the peak of -"^^CO. These different characteristics of 
^^CO and H I can be explained within the frame of TASS model where ^^CO clouds mostly 
correspond to the high-density region just behind the spiral shock. We have compared the 
I — V distribution of ^^CO around 6 ~ 0° with the theoretical prediction (shock ridge) and 
found them to be well matched. Hence, it is likely that the distances to the Tycho SNR and 
to the molecular clouds surrounding it are comparable. 

3.2. Non-radiative Ha filaments 

Optical filaments have been found in Tycho (Smith et al. 1991). The filaments are 
located along the NE boundary of the remnant, where interaction with a molecular cloud is 
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likely (Figure 3). The emission is predominantly hydrogen B aimer emission and is believed 
to originate from a coUisionless nonradiative shock, characterized by broad and narrow emis- 
sion components whose velocity distribution represent that of postshock and preshock gas, 
respectively (Chevalier & Raymond 1978; Chevalier et al. 1980). The expansion velocities of 
these filaments, which are compatible with radio expansion rate at this region (~ 0.15"/yr, 
Reynoso et al. 1997), is relatively small (~ 30 %) compared to the other parts of the rem- 
nant. The implied large preshock density suggests its association with the molecular cloud. 
Ghavamian et al. (2000) obtained high-resolution echelle spectrum of knot g (Figure 3) and 
their observed central velocity is i'lsr = —53.9 ±1.3 km s~^. For comparison, most of the 
^^CO emission around Tycho is between —67 and —60 km s~^, and the ^^CO line profile near 
knot g, fitted with Gaussian, gives a line center of -62.50 ± 0.01 km s"^ and a FWHM of 
2.27±0.03 km s^^. Therefore, there is a non-negligible velocity difference of 8.6± 1.3 km s^^ 
between the H alpha filament and the ^^CO cloud. Even if we consider the FWHM of ^^CO 
as a possible error, the difference is greater than 6 km s^^ which is significant. The velocity 
of the narrow Ha component of a non-radiative shock should represent that of the gas enter- 
ing the shock. This imphes two possibihties. One is that the gases responsible for ^^GO and 
Ha have initially different velocities. This could be due to a turbulent motion behind the 
spiral shock or some kinematical effects associated with the evolution of the progenitor star 
or the SN explosion itself. But these seem rather unlikely because there is virtually no ^^GO 
emission above —60 km s~^ and the actual velocity difference considering a projection could 
be much larger (see next section). The remaining possibility is that the gas experiences an 
acceleration before the excitation by the shock. We further discuss about it in the following. 

A handful of remnants show non-radiative shock emission like Tycho (Smith et al. 1991, 
1994; Ghavamian et al. 2001; Sollerman et al. 2003). The width of the narrow line has been a 
puzzle since it is unusually large for neutral hydrogen and does not vary much from remnant 
to remnant. For most of the remnants, it ranges between 30 and 50 km s~^, while SN 1006, 
which is the fastest (~ 3000 km s^^), has a slightly smaller value of about 20 km s^^. Various 
physical mechanisms have been proposed to explain this line width (Smith et al. 1994). 
Although the GR precursor and/or fast neutral precursor scenario has been regarded as the 
most likely case (Smith et al. 1994; Ghavamian et al. 2000), firm evidence is yet to be found. 
Some of the proposed mechanisms, which include a interaction with heavy particles (like 
protons), also predict a systematic acceleration of the gas before their excitation. We propose 
that the velocity difference between Ha and ^^CO can be explained by this gas acceleration. 
Observed peak velocity of the broad Ha component is red-shifted by 132 ± 35 km s~^ from 
that of the narrow one, indicating that the shock is not completely edge on (Ghavamian 
et al. 2001). With a shock velocity Vg — 2000 km s~^ (fluid velocity ~ 1^^), viewing angle 
is estimated to be about ^ ~ 5°. If the gas is accelerated, the Ha narrow component will 
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be red-shifted from the unperturbed medium, which is quahtatively consistent with our 
interpretation where ^^CO typifies the velocity of the unperturbed medium. The observed 
velocity difference of 8.6 ± 1.3 km in wlsr then implies that the actual acceleration 
is about Av = 98 ± 30 {vg / 2000 km s~^) km s~^. One possibility to explain the velocity 
difference is a momentum transfer by collisions with fast protons after crossing the shock. 
But, as Smith et al. (1994) has showed, a lifetime of the hydrogen atom is very short and 
only a moderate amount of acceleration (< 10 km s"^) is possible. We consider two preshock 
mechanisms which might account for this acceleration. 

Broadening of the narrow component is well explained by a CR precursor (Smith ct al. 
1994) and acceleration of the preshock gas is also theoretically well expected (see Blandford 
& Eichler 1987). For a CR dominated shock, it is suggested that the preshock gas can 
be accelerated to the significant fraction of the shock velocity. Boulares & Cox (1988) has 
applied this shock model to the B aimer-dominated filaments in Cygnus loop. According to 
their model, the preshock gas is accelerated to velocities close to the shock velocity. For the 
case of Tycho, the smaller amount of acceleration, compared to the shock velocity, indicates 
that the cosmic ray contribution is not dominant. This is consistent with recent X-ray 
results (Hwang et al. 2002) where no significant temperature drop due to a possible cosmic 
ray acceleration is observed. The radio spectral index (5*^ oc t/~") of Tycho is known to be 
steeper (a ~ 0.6 at u > 1 GHz) than what is expected from a test particle theory (0.5), 
and it is suggested that Tycho may have somewhat strongly modified shock (Reynolds & 
EUison 1992; Volk et al. 2002). But, the spectral index varies over the remnant (Katz-Stone 
et al. 2000) and knot g region (filament I in Katz-Stone et al.) has the flattest spectral index 
of 0.44 ± 0.02, which might be another indication for the interaction with molecular clouds 
(Koo ct al. 2001). Although Tycho may have strongly modified shock in a large scale, it 
does not seem to apply to areas interacting with molecular clouds. 

There has been little theoretical study on the fast neutral precursor compared to the 
cosmic ray precursor. Lim & Raga (1996) calculated the Ha profile from a 225 km s~^ 
non-radiative shock with a fast neutral precursor, but due to their narrow parameter range, 
direct comparison with our observation is not available and there is no indication of the 
velocity shift of the narrow component in their line profiles. Calculations for higher shock 
velocity and lower pre-ionization fraction is necessary. 

4. CONCLUSION 

The study of Balmer-dominated filaments in SNRs has given us much information about 
the shock structure. The velocity of an unperturbed preshock medium is an invaluable 
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observational ingredient for fully understanding underlying shock structure. In Tycho, a 
velocity shift between an ambient molecular cloud, which is believed to be the unperturbed 
preshock gas, and the Ha filaments has been observed. Although the preshock acceleration 
was searched for previously, no evidence has obtained before our observation (e.g. Ellison 
et al. 1994). Our observation serves as the first direct evidence for the acceleration of preshock 
gas by a cosmic ray (or fast neutral) precursor. We, however, note that our conclusion is 
based on the single velocity measurement and more observations are needed to test it. 

There has been not much study on the structure of collisionless shocks propagating into 
molecular clouds. But most of the vohimc of a molecular cloud may be occupied by low 
density (< 10 cm^'') intcrclump medium, and the structure would not much differ from that 
in atomic medium. The molecules are supposed to be dissociated within a precursor either 
by radiation from postshock region or hot electrons. The dissociation rate of H2 is only 
slightly larger than ionization rate of the hydrogen atom (Yan et al. 1998; Le Teuff et al. 
2000), and complete dissociation of H2 may not happen since they occupy more dense region. 
The velocity of the transmitted shock into a molecular clump can be sufficiently low to be 
radiative, but currently there is no sign of radiative shock in Tycho, presumably because 
the shock is propagating into the low-density outskirts of the molecular cloud. Future study 
should address these issues together with the nature of cosmic ray acceleration for a high- 
velocity (> 1000 km s~^) shock propagating into a molecular cloud. 
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Fig. 1. — Integrated channel map of ^^CO from FCRAO data (—68 < v < —59 km ). 
White contour for 1420 MHz radio continuum. Black boundary represents Nobeyama ob- 
served region. 



Fig. 2. — ^^CO channel map of Nobeyama observation with Av = 1 km . Contour 
represents 1420 MHz radio continuum. The last panel shows average spectrum of total 
observed region. The peak brightness temperature of the average spectrum is 0.83 K. 



Fig. 3. — Pseudo color image of Integrated channel map (—63 < v < —60 km s ^). White 
contour for 1420 MHz radio continuum, red for Ha (Ghavamian et al. 2000), and blue for 
hard X-ray (4 keV < hu < 6 keV). 
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